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Abstract 



^^ ' We construct a simple, robust model of the chemical evolution of 

^^ . galaxies from high to low redshift, and apply it to published observa- 

f^ ' tions of damped Lyman-alpha quasar absorption line systems (DLAs) . 

f^ , The elementary model assumes quiescent star formation and isolated 

^D ' galaxies (no interactions, mergers or gas flows). We consider the influ- 

r~| . ence of dust and chemical gradients in the galaxies, and hence explore 

pi i' the selection effects in quasar surveys. We fit individual DLA systems 

to predict some observable properties of the absorbing galaxies, and 
^ ' also indicate the expected redshift behaviour of chemical element ratios 

C/3 , involving nucleosynthetic time delays. 

^ ■ Despite its simplicity, our 'monolithic collapse' model gives a good 

account of the distribution and evolution of the metallicity and column 

density of DLAs, and of the evolution of the global star formation rate 

and gas density below redshifts z '^ 3. However, from the comparison 

C^ ' of DLA observations with our model, it is clear that star formation 

rates at higher redshifts (z > 3) are enhanced. Galaxy interactions 
and mergers, and gas flows very probably play a major role. 

1 Introduction 

Observational studies of the formation and evolution of galaxies are making 
great strides forward. Unfortunately, understanding of the wide range of 
physical processes involved in this key astrophysical problem is far from 
complete. There is a pressing need for a simple yet robust framework in 
which to interpret the new observations. 

The prevailing view is that galaxies form in dark matter halos that trace 
primordial density fluctuations, but subsequent mergers, gas flows and star 
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formation certainly play an important role in galaxy evolution. Great efforts 
have been made in the 'semi-analytic' approach to modelling galaxy evolu- 
tion, as summarised by ( Baugh et al., 199^ ). This sophisticated technique 



uses statistical methods to follow the growth of dark matter halos, and phys- 
ically motivated rules to describe the gaseous and stellar processes within 
them. Ah initio calculations of galaxy evolution in hierarchical cosmolog- 
ical models are constrained by local galaxy luminosity functions, enabling 
statistical and cosmological predictions such as galaxy survey redshift dis- 
tributions, and the cosmic star formation history, to be made. 

Our approach is different. We have produced a simple, robust model 
for galaxy formation and evolution intended to meet the constraints of the 
observed chemical evolution of the Universe. In particular, in this paper 
we investigate the detailed gas phase abundances of quasar damped Lyman 
a absorption line systems (DLAs) to z ~ 4-5 (defined as those with neu- 
tral hydrogen column densities above log (A'^^/ cm~^) ~ 20). We improve 
upon elementary global calculations (e.g. (|Pei and Fall, 1995t pi^dmunds and 



Phillipps, 1997] )) by modelling the enrichment of gas in separate gravita- 
tional potential wells. We adopt as far as possible the 'best' descriptions of 
the different physical processes involved, based upon reliable observations 
and theory. Our eventual aim is to predict the outcome of future observa- 
tions - particularly in the IR and submm, with very large telescopes and 
adaptive optics - that will allow chemical analysis of directly detected indi- 
vidual galaxies or galaxy fragments at high redshifts. 

The DLA systems have been linked to proto-galactic clouds ( ([Wolfe 



et al., 19861) ), low surface brightness galaxies (( [Jimenez et al., 1999| )), the 
outer regions of galactic disks (( [Bergeron, 1986| ; [Bergeron and Boisse, 1991 



Bteidel, 1995| )), sub- galactic star forming regions, and dwarf galaxies (( jYanny 



and York, 1992D ). DLAs are currently known at redshifts 0.1 < z < 4, which 
probably spans the era of major galaxy formation, and have observed HI col- 
umn densities up to log (Nhicui'^) ~ 21.7. We do not attempt to address 
the precise nature of the absorbers, but treat them as clouds of gas (with 
associated stars and dust) irrespective of what kind of (proto)-galaxy they 
may represent. We have compiled a database of redshifts, HI column densi- 
ties and abundances of all the DLA systems with published measurements. 
We have produced synthetic surveys for overall statistical comparisons, and 
models of individual systems, and hence predicted directly observable sys- 
tem properties. 

The layout of this paper is as follows. In § |2|, we outline our initial 
model for the formation and evolution of the galaxy systems, using the 
Simple chemical evolution model. In § |^, we discuss the dust optical depth, 
and add a simple representation of galaxy internal structure to our model, 
in the form of exponential disk gradients in stars, gas, metallicity and dust. 
In § ^, we consider the opacity of galaxies and possible effects upon the 
quasar optical surveys. In § 0, we use the Holmberg relation for galaxy 
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luminosity and size to make a preliminary calculation of the evolution of the 
galaxy population as a whole. In § |6|, we compare detailed model predictions 
with observations of DLAs, including the chemical evolution of individual 
elements, and study other relevant data such as the properties of Local 
Group dwarf galaxies, and the star formation history of the Universe. In 
§ 1^, we take a critical look at the limitations of our initial model, and how 
these might affect our conclusions. In § |^ we summarise our conclusions. 
(Unless otherwise stated, our calculations assume Hq = 100/ikms~^ Mpc~^ 
with h = 0.7, and Vt = 1.) 

2 A Galaxy Chemical Evolution Model 

We wish to focus upon chemical evolution, and so we construct a simple, 
yet physically reasonable model of a galaxy, without imposing particular 
structural or dynamical properties. We consider a spherically symmetric 
cloud of gas which, at a redshift determined by its density, breaks away 
from the Hubble flow and collapses to form a gravitationally bound system. 
Given a redshift, the density required for a region to collapse is (( [Peebles! 
199^ )) 



0.3(z + 1)3^2 Om(<5 + 1) ,.,, 

P' = 2G ^^^ 

where the density contrast over the mean background density 5 = 5.6, 
and Q.M is the mass density of the Universe (for cosmologies with Q\ = 0) . 
Guided by detailed 3D numerical self-gravitating hydrodynamical calcula- 
tions of the formation of galactic systems {e.g. ( phurches, 1999| )), we assume 
that the cloud collapses by a factor of 10 from an initial radius Rq over a free- 
fall time, halting at i? = 0.1i?o and evolving at constant radius thereafter. 
The free-fall time of the cloud, tg is 

where pQ is the density of the cloud as it breaks away from the Hubble 
flow at the turnaround redshift, ztum- The radius of the cloud, R, is given 
(for t < t^ ) by 

^-^(^-ijj- (3) 

We model clouds with total gravitational mass in the range 10^ < 
M/Mq < 10^2. We fix the baryonic fraction to be universally 10%,, so 
that each cloud contains 90% dark, non-baryonic mass which does not form 
stars, nor contribute to the gas column density. For clouds of uniform den- 
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sity, the initial cloud radius (at the onset of collapse) is therefore a known 
function of cloud mass and the turnaround redshift (Eqn |l|) . 

We use a Schmidt star formation law (( ^chmidt, 1959| )), with the star 
formation rate a power-law function of Pgas-, the gas volume density, given 

by 

where V is the volume star formation density, and k and q are constants 
which determine the star formation efficiency. Kennicutt has shown that 
such prescriptions are a good representation of star formation over a wide 
range of physical conditions in present-day galaxies (( [Kennicutt, 1988| ; Ken- 



nicutt, 199^ )). We fix the slope at a convenient value, q = 1.5, guided 
by Kennicutt's prefered slope 1.4 it 0.1. We fix the absolute efficiency 
K = 2.5Mq ~^'^pc^'^Gyr~^, which gives reasonable present-day gas/star 
mass ratios across a range of initial model galaxy cloud masses (confirmed by 
detailed self-consistent numerical galaxy models ((Churches, 1999|))). The 



star formation rate in all our calculations therefore varies only as a function 
of the gas density. 

The total baryonic mass of the model galaxy cloud is M5 = O.IM. Ini- 
tially, just after turnaround, the cloud baryonic mass is purely gas, so the 
mass of gas in the cloud Mg = M^, and the mass of stars, Ms = 0. Star 
formation proceeds during collapse and subsequent galaxy evolution and we 
assume that a fraction, a, of the gas which is consumed by star formation re- 
mains in long-lived stars and stellar remnants. Therefore, the rate of change 
of the mass in stars is given by 

^ = '-a.R^T (5) 

dt 3 ^ ' 

at all times. Substituting for the volume star formation rate, F, from 

Eqn ^, and the current radius of the cloud, i?, from Eqn ^ and noting that 

Mg = Mb- Ms we find 



^"^■^ -^.nil-tffr^^Mb-Msf' (6) 

^0 



dt 
for the cloud collapse phase, and 



dMs 1000a , ,3/2 

— - = — ^ K. (Mu - uyi^ 

dt y 



jj^ K (Mb - MsY'' (7) 



for the constant radius cloud evolution phase, where Vq is the initial 
volume of the model galaxy cloud. Both these equations may be solved 
explicitly to give 
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Mg 



\ -2 

2aK 



yo 



"T^{tff-tffyHtff-tf')+M^^/'] (8) 



and 




^9=4^ (0-063 td + 2% ) + M-^'^ ) (9) 

for the collapse and constant radius phases respectively, where td is the 
time since the end of the collapse phase. 

2.1 Chemical Evolution: Testing the Model 

The fundamental aim of this work is to make testable predictions about the 
chemical evolution of galaxies at high redshift, and in the first instance, to 
understand the properties of high column density (damped Lyman alpha) 
quasar absorption line systems. Therefore, we need to calculate the basic 
observable properties of galaxies which appear as DLAs. These are: the ab- 
sorption redshift Zabs, the gas phase metallicity, Z (mass fraction of metals in 
the interstellar medium) and the column density in neutral hydrogen, Nhi , 
in the ISM. The absorption redshift is simply determined by the turnaround 
redshift of the model galaxy cloud, its age, and the cosmology. We can also 
readily calculate an average metallicity and column density for any cloud. 

The metallicity can be estimated from the Simple Model of 'closed box' 
chemical evolution (( pagel, 1997] )). The standard Simple result is that the 
metallicity, Z, follows the relation 

Z = -plnf (10) 

where p is the yield (the mass of heavy elements in the processed gas 
returning from stars to the ISM per mass locked up in long-lived stars and 
stellar remnants), and / = Mg/M is the fraction of baryons in gas (hereafter, 
the 'gas fraction'). Throughout our calculations, we assume a constant, 
uniform heavy element yield, p. 

The column density of neutral hydrogen, Nhi can then be estimated 
from the cloud gas mass Mg (allowing for primordial helium, and assuming 
a typical cloud size R) as Nhi ~ Mg R/{^mH){^Tr R^), giving 

16mH vr R^ 

where nin is the mass of a hydrogen atom. This is an upper limit, since 
molecular hydrogen will form in the densest regions, and there will be some 
ionisation. 

We now have a basic framework of galaxy formation and evolution, and 
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in particular, a robust formalism to study the chemical enrichment history 
of any given model galaxy cloud. We next need to study the dust optical 
depths that can be expected (§ ^) (which may affect quasar surveys (§ [4.1D ), 
and to allow for internal structure within our galaxy clouds (§ |3.2| ). 

3 Gradients and Dust 

A significant fraction of the heavy elements in our Galaxy are locked up in 
the solid phase, particularly as small dust grains in the ISM. We can predict 
the evolution in the opacity of our model clouds due to absorption by these 
grains, Td, in terms of Nhi and gas fraction, /. If a cloud is to be detected 
as a DLA, we must first detect the background quasar, often in a survey at 
blue wavelengths; the dust must not redden and obscure the quasar beyond 
survey and spectroscopic magnitude limits (( |Pei and Fall, 1995 )). 



3.1 Average Dust Opacity 

The optical depth of a cloud due to dust r^ = N^ 0"^ where a^ ~ vrr^ is the 
typical grain absorption cross-section (r^ is a typical grain radius) and iV^ 
is the column density of dust grains. For a typical line of sight, the column 
density of dust N^ ~ n^ R, where n^ is the number density of grains and R 
is the characteristic size of the cloud. Therefore, the dust optical depth can 
be estimated from the assumed typical grain properties, and global cloud 
properties 

Td = ndRTrrj. (12) 

If Md is the total mass of dust in the cloud then we can write (( [Edmunds| 
and Eales, 19981) , using our symbols) Md = rjZMg, where Z is the mass 



fraction of metals in the gas phase, rj is the fraction of metals which have 
condensed into dust and Mg is the cloud gas mass. We can also write Md = 
nd rrid V, where V = ^tt R^ is the volume of the cloud; if the typical density 
of a dust grain is pd then the typical mass of a dust grain is nid = ^irr^pd- 
Therefore, equating our two expressions for Md, 

vZMg = nd{^7TfR^rjpd. (13) 



Thus, combining Eqn 13 with the expression for dust opacity (Eqn p^) , 
and using the relationship between gas fraction and metal mass fraction for 
the Simple Model (Eqn p!o| ) , the optical depth due to dust may be written 
as 

9Mo?7nln4 
r, = ,/', ^ (14) 

iQirR'^Vd Pd 
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The cloud gas mass and column density are related by Eqn 11 so we can 
eliminate cloud gas mass Mg in the expression for optical depth (Eqn |lj) . 
The typical dust optical depth in a model galaxy cloud is therefore 

Td = ^, (15) 

rdPd 

Theoretically, we calculated the gas fraction, /, and hence, the column 
density, Nhi ■ Observationally, we measure column density and the metal- 
licity. We can therefore use Eqn ^ to infer average dust opacities for model 
galaxy clouds anc? observed systems, if we make some assumptions about the 
properties of the dust. Our knowledge of the Galaxy suggests that taking 
as representative values ry = 0.5 (the fraction of metals in dust), p ~ 0.01 
(the local stellar yield, Q.^Zq ), r^ ~ 0.1//m (the radius of a typical 'classical' 
dust grain) and pd ~ 2000kg m~'^ = 2gcm~^. With these values, the typical 
dust optical depth of a cloud is 

Trf ~ 0.4iVi^/,2i In - (16) 

where Nuj = Nhi^21 x lO^^cm"^. This prescription indicates that a 
model galaxy cloud representative of our Galaxy, with a gas fraction / ~ 0.1 
and an HI column density given by Nhi,21 ~ 1 will have a average dust 
optical depth, r^ ~ 1. It is possible that the fraction of heavy elements 
locked up in dust in the ISM (77) is a function of galaxy evolution, but initial 
estimates suggest that this is not a significant effect for purposes of the 
present calculation (( [Edmunds, 2001 )). 



3.2 Gradients 

Our calculations so far have implicitly assumed that the dust is uniformly 
distributed. We now allow for radial exponential gas, metallicity and opacity 
gradients, as found in disk galaxies. 

We assume that the metallicity distribution does not vary azimuthally, 
and that Z{r) can be parameterised as log Z{r) = a — br, where a is the cen- 
tral metallicity of the disk and the slope, b ~ 0.2 dex/scalelength (( parnefrt 
199^)). Equivalents, 



Z = Zo e-°-^^"/^* , (17) 

where h^, is the stellar exponential scalelength of the disk. We calibrate 
using the Galactic measurement Zjsm = 0.017 at the Solar circle, r ~ 8kpc. 
Then, if Mb,mw is the absolute blue magnitude of the Galaxy, the central 
metallicity of a galaxy with absolute blue magnitude Mb, is 

Zo = o.027e-°-^^(^^s-^^s-*^^). (18) 
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For a typical L* galaxy, the neutral hydrogen is also distributed in an ex- 
ponential disk, and the HI scalelength and the stellar scalelength are roughly 
the same: h ~ Hhi ~ /i* ~ 7 kpc. Then Nniir) = Nhiq e~^'^ and the opac- 
ity gradient is therefore (from Eqn |l^) 

, , 1.3 X rimij , . , . , . 

r(r) = l_JLNHi{r)Z{r). (19) 

Pd^d 

The central HI column density, and opacity, are uncertain because con- 
version to molecular hydrogen may promote the formation of dust, especially 



in cooler regions ((Valentijn and van der Werf, 1998|)). Substituting in nu- 



merical values as in the previous section and normalising to the column 
density at one scalelength (meaning that Nhi^^ = 4:NHi^-^{h) x 10^^ cm~^) 
gives 

T{r)r. 1.5 NHi,Ah)e-'-''^'/'' (20) 

where the column density measure, NHi^h,2i is directly calculated in our 
model. This relation predicts an average optical depth at the Solar circle in 
our Galaxy of r^ ~ 1, which is again consistent with azimuthal averages in 
external galaxies ( ([Domingue et al., 1999| ), (|Kuchinski et al., 1998]) ). 

4 Line-of-sight Optical Depth 

The analysis of galaxy opacity gradients presented above suggests that, for 
a present-day L* galaxy, the disk optical depth becomes low enough to allow 
the detection of an underlying L*^ quasar in typical optical surveys at ~ 2 
disk scalelengths; hence ~ 50% of the disk surface area will be 'optically 



thick' to -L^ quasars (see § iA). We can therefore set an upper limit on 
the fraction of the sky which is covered by the optically thin regions of disk 
galaxies. The probability P that a random line of sight passes through a 
given region of a galaxy varies with redshift as 

dP c(l + zf 

where n is the number density of galaxies, and a is the geometric cross- 
section of the relevant region of the galaxy (( Peebles, 1993| )). In general 



n 



and a will be (probably decreasing) functions of redshift, but to determine 
a firm upper limit, we can assume that n and a are constant in redshift. 
In an Einstein-de Sitter universe, E{z) = (1 -|- z)^'^, and so the probability 
that a random sightline intersects the outer regions of a galaxy where there 
is insufficient dust to obscure an L* quasar, yet sufficient neutral hydrogen 
to give rise to a DLA is 
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P='^((l + z)"'-l)^ (22) 

We take n ~ 0.008 Mpc"'^ ((jBinney and Merryfield, 1998|)). From our 



opacity gradient calculations (§ 3^), the optically thin cross-sectional area 
of a typical face-on disk is ~ 600 kpc^. Therefore, our firm upper limit on 
the probability of a line of sight intersecting such a region is ~ 10%, for 
a quasar at an emission redshift Zem ~ 3. Of all the quasars that have 
been discovered, ~ 10% of quasars which have been spectroscopically sur- 
veyed have damped absorption features (( [Wolfe et al., 199^ )). This very 
simple analysis would therefore suggest that all DLAs are caused by the 
outer regions of L* galaxies. How is this conclusion affected by our a priori 
assumptions? We have assumed that all galaxies are oriented face-on, which 
gives an upper bound to the optically thin cross-section. Implicit in the 
above reasoning is the assumption that the population of quasars forms a 
shell at high redshift {e.g.z = 4), maximising the intersection probability. 
We assume no evolution in the number density of the disk population, but if 
the epoch of disk formation falls at z < 4, the number density at earlier times 
will be lower. All these factors lead to overestimation of the probability of 
intersection. Hence, the covering fraction on the sky of suitable regions of 
L* galaxies (optically 'thin', 'high' Nhi ) is probably insufficient to explain 
the frequency of DLAs. The DLA progenitors are probably a heterogeneous 
population, with some absorption features forming in the outer regions of 
L* galaxies, but others forming, for example, in the central regions of dwarf 
galaxies, or in low surface brightness galaxies. 

4.1 Quasar Selection Biases 

Intervening galaxies with the dust properties discussed above (§ |3|) may cause 
selection effects in quasar and DLA surveys. These effects will depend upon 
quasar survey wavelengths and magnitude limits, and upon spectroscopic 
follow-up magnitude limits, and upon the density and luminosity distribu- 
tion of quasars. 

The form and evolution of the luminosity function of optically selected 
quasars is quite well known. The shape of the quasar optical luminosity 
function (QOLF) ^ (L, z) can be well parameterised as a double power-law 
(( Boyle et al., 1988[ ); see Fig[l|a). The optical QLF has just one character- 



istic feature, the 'knee' ^*,M^, which enables us to distinguish between 
luminosity and density evolution in the quasar population. The QOLF 
moves monotonically towards higher luminosity as a power law in redshift, 
oc (1 + z)^ , where the evolution parameter is A; ~ 3.4 (see recent deter- 
minations by ([Franca and Cristiani, 1997 ; Maloney and Petrosian, 1999[ ), 



and Fig [l[b). Interestingly, this intrinsic brightening of the characteristic L* 
quasar rather closely balances the effects of cosmological dimming (Fig ||c) . 
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[INSERT FIGURE 1 HERE] 

The QOLF is calculated from quasars discovered in magnitude-limited 



blue surveys {e.g. the LBQS, ( Hewett et al., 1995 )). Quasar numbers are 



dominated by those fainter than the characteristic luminosity L* at all red- 
shifts. The level of foreground obscuration which will dim an L* quasar 
below the typical survey magnitude limit is subject to uncertainties in the 
adoption of cosmological and quasar evolution models. However, for median 
values of all the parameters, in the redshift range z > 0.8 it would take fore- 
ground extinction of only r ~ 0.5 to push an L* quasar below an apparent 
magnitude limit msi^^ = 19 (Fig ||d). 

In Fig. y, we compare the typical opacities calculated from metallicity 



and HI column density in § 3.1 with observed properties of DLAs and nearby 
galaxies. The majority of known DLAs have r < 0.1, which will not cause a 
significant number of quasars to be lost from optical quasar surveys. How- 
ever, it is intriguing to note that the DLA data points cluster along a line of 
constant opacity, approximated t^ = 0.1. Despite other existing evidence to 



the contrary {e.g. Ellison et al. (in preparation), and see § 6.2.1, the simplest 
interpretation of this analysis is that many DLAs do exist with higher opac- 
ities, but are missing from our catalogues because they obscure background 
quasars beyond optical quasar survey magnitude limits. 

[INSERT FIGURE 2 HERE] 



5 The Evolution of the Galaxy Mass Distribution 

To make statistical statements about the DLA population, we need to under- 
stand the form and evolution of the mass function of galaxies. In our model, 
the metallicity evolution of a galaxy is only a function of its turnaround 
redshift, which uniquely determines the gas density (Eqn ||). To break 
this degeneracy, we need a relationship between mass (or equivalently, ra- 
dius), and turnaround redshift. We can obtain such a relationship (following 
( Phillipps and Edmunds, 19961) ), from the Holmberg radius- luminosity rela- 



tion (( [Holmberg,"l97^ )): 



Rhb {L) « L^ (23) 



where Rhb is the Holmberg radius, L is the luminosity, and the parame- 
ter s ~ 0.4. Assuming a constant mass-to-light ratio, the Holmberg relation 
implies a mass-radius relation Rq = AM^ {A is a constant, Rq is the initial 
cloud radius) and hence a mass-density relation. We can therefore calcu- 
late the (unique) turnaround redshift for each galaxy mass. Since we are 
not including galaxy mergers in this calculation, and we know the present 
day galaxy luminosity function, we can construct a history of galaxy forma- 
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tion in the universe. Large numbers of low mass dwarf galaxies form first, 
with higher mass galaxies forming at progressively lower redshifts in ever 
decreasing numbers. 

We can substitute our mass-radius relation into the expression for the 
turnaround redshift (Eqn [l|) , and re-arrange to obtain mass as a function of 
turnaround redshift: 

The overall present-day space density of galaxies of mass M in mass 
interval dM is well modelled by the Schechter luminosity function 

cI>(M)dM = CD* (^_ j exp(^-— jdM (25) 

where $* and M* are the characteristic galaxy number density and mass, 
respectively. The amount of mass per co-moving unit volume bound up in 
galaxies which form at a given redshift is given by 



KA ( \ ^*M* II {zturn + lfA^Hla{5 + l 
Mgal[Zturn) = ^ M exp 1 - 




^ 2G 

(26) 
Unfortunately, this result is extremely sensitive to the slope, s, chosen 



for the Holmberg relation (Eqn 23). A value of s = 0.4 produces a uni- 
verse populated by galaxies that are too young and metal poor, whereas 
s = 0.42 (Holmberg's own value) produces a universe populated by over- 
evolved galaxies. The necessity for fine tuning of this parameter is unsat- 
isfactory. Another serious problem is that Holmberg's original tight corre- 



lation between radius and luminosity (( Holmberg, 1975 )) is much less con- 
vincing in large modern datasets (Gavazzi, private communication). There 
is considerable scatter, and selection biases which are caused by the diffi- 
culty in discriminating between stars and high luminosity compact galaxies 
on one hand, and insensitivity to low surface brightness galaxies on the 
other, suggesting that the Holmberg relation may be no more than a selec- 
tion effect. We may be able to side-step these problems entirely by using 
the Extended Press-Schechter theory in our galaxy formation models, and 
including hierarchical galaxy evolution through mergers (see § |3). We will 
consider this model in our next paper (Baker et al. 2001). Nevertheless, 
our simpler "Holmberg" model is surprisingly successful in explaining the 
observed chemical evolution of DLAs, as we now discuss. 
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6 Observations Confront Theory 

We are now ready to confront our preliminary theoretical predictions with 
the observed properties of damped Lya absorber systems from the literature 
{e.g. absorption line survey data ( Rao and Turnshek, 2000| ; Storrie-Lombard 



et al., 1996a; Lanzetta et al., 1995| |Junkkinarien et al., 1991 ; Lanzetta et al.. 



1991 



Wolfe et al., 1986| ), reports of metallicity measurements ( |Centurion| 



et al, 2000|; pcttini et al., 2000|; [Boisse et al., 199^; [Centurion et al., 1998|; 



Molaro et al., 1998; Matteucci et al., 1997; Pettini et al., 1997; Lu et al., 



19961 ; pmette et al., 199^ ; jLu et al., 1995|; [Meyer et al., 1995| ; pettini et aL 



;pu 



1995 ; Pettini et al., 199^ ; [Meyer and York, 1992 ), and metallicity measure- 



ments with 8m-class telescopes ( putram et al., 1999[ ; Pettini et al., 1999a 
Prochaska and Wolfe, 1999| ; [Prochaska and Wolfe, 1997 ; Prochaska and 
Wolfe, 1996| )) The free parameters of the model clouds are the turnaround 



or formation redshift (which fixes the density) , and the mass of the collaps- 
ing halo, which together uniquely determine the initial cloud radius. After 
commenting upon metallicity indicators (§ 3.1), we will now discuss our 
simulations of DLA survey data (§ [6.2[) in some detail, including (§ [6.2.2 ) 
implied global properties (star formation, gas density). In § 3.3, we examine 
fits to individual DLAs, then (§ [6.4[ ) briefly compare the properties of DLAs 
with dwarf galaxies in the Local Group. In § |6.5[ then consider the chemical 
evolution of time-delayed species. 



6.1 Gas Phase Metallicity Indicators 

Zinc has been considered as a fair and useful indicator of the gas phase 
metallicity because Zn is measured to follow Fe in Galactic stars, suggesting 
that Zn might track the stellar nucleosynthesis of Fe in Type la supernovae 
(( Arnett and Thielemann, 1985 ; Thielemann and Arnett, 1985[ )). However, 
recent measurements in metal poor and thick disk stars find a ratio of [Zn/Fe] 
~ 0.1 ( ([et al., 2000[) ), perhaps consistent with an origin in Type II super- 
novae (( Arnett, 1996D ). Whatever its stellar origin, Zn is expected to be 



far less depleted onto dust grains than Fe (( Pettini et al., 1999b )). Unfortu- 
nately, the weakness of the observationally accessible Znll transitions, whilst 
circumventing problems of saturation, means that Zn is extremely difficult 
to measure at Zahs > 3 (( [Prochaska and Wolfe, 2000| )). Current measure- 
ments of DLA abundances show < [Zn/Fe] >~ 0.6, whereas <[Zn/Cr]>~ 0.3, 
implying depletion factors of 4 and 2 for Fe and Cr respectively if Zn is un- 
depleted. As accurate measurements of [S/Si] become available for DLAs 
with [Zn/Cr] < 0.3, it will be possible to test for metallicity-dependent de- 
pletions (( Pettini et al., 2000| )). In the meantime, we too will use Zn as our 
primary metallicity indicator. 
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6.2 Simulated DLA Surveys 

The principle observations which we are attempting to explain are the red- 
shifts, column densities and metallicities of DLA systems. To compare our 
model with these observations, we created artificial surveys of 1000 DLAs, 
including overall observational selection effects in redshift, column density 
and metallicity. We consider our model in two simple yet representative 
cosmologies: the Einstein-de Sitter (EdS; ^m = 1)) and a low density cos- 
mology with J7a/ = 0.3, r^fc = 0.7; JIa = and h = 0.7 in both casesQ We 
have sampled the model DLA galaxies at random impact parameters along 
the gas and metallicity radial gradients. The results, in terms of Zabsi^Hi 
and [Z/H], are summarised in Figures § and g, with the observed data points, 
and simulated survey number density contours. Absorption redshift is the 
random input parameter for each model DLA system, and we use the red- 
shift selection function shown in Fig. |3|, representing the combined redshift 
survey path lengths X{z) of the DLA surveys of ( [Lanzetta et al., 1995 ; 



Btorrie-Lombardi et al., 1996b ) and ( Rao and Turnshek, 2000 )) 



It is clear that, in the EdS cosmology, our model is a poor description 
of the data for redshifts z > 1.8 (see Fig. |^, and we will concentrate on our 
other representative model, the 'low density' cosmology, in the rest of this 
discussion. 

[INSERT FIGURE 3 HERE] [INSERT FIGURE 4 HERE] 

For HI column density, the observational detection threshold is set at 
log {Nhi cm^^) > 20.3 in all statistically complete DLA surveys. For metal- 
licity, the sensitivity limit depends upon HI column density, and there is no 
uniform Zn equivalent width limit applied to observational surveys. There- 
fore, we have estimated the effective metallicity selection function from the 
surveys of Pettini et al. to be [Zn/H]mm ~ —2.5 log [NHicm.'"^) +50 (( pet- 



tini et al., 2000| ) and references therein; ( prochaska and Wolfe, 2000| ) quote 
a similar combined limit, that log (A^/f/ cm~^) + [Zn/H] > 19.0). These 
selection functions are indicated in Fig. ^ as dashed lines, and have been 
folded into Figs. ^ and |6[ 

6.2.1 Cumulative distributions 

The cumulative distributions of redshift, column density and metallicity for 
the low density cosmological model are shown in Fig. 0. The smoother line is 
the distribution of simulated DLA systems and both curves are normalised 
to show the fraction of 'surveyed' systems included up to a given parameter 
value. The observed and model cumulative redshift distributions agree well. 



^We assume that the contributions of matter, cosmological constant and curvature are 

such that ^M + riA + rife = Q.total = 1- 
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The column density measurements are subject to a Malmquist type bias 
but for the column density distribution function calculated by ( ^torrie- 



Lombardi et al., 1996a) and the existing absorption redshift survey path 
lengths, X{z), only about one spurious extra DLA system is expected to be 
scattered into the lowest column density bin. Again, the model and observed 
column density distributions agree well. 

The cumulative metallicity density distributions (Fig. ^c)), show that we 
have a deficit of low metallicity systems ([Zn/H] < —1.3). This may relate 
to star formation at high redshifts, a point which we return to in § |6.2.2 . 



The predictions in the low density cosmology are clearly a significantly 
better representation of the data that those in the EdS cosmology, especially 
at redshifts z < 3 where the Zn measurements exist. The well-known ob- 
served deficit in high HI column density systems {log [N h i cvot'^^ > 21.0) 
at high redshifts (z > 3.5) is reproduced (Fig. ^), because the sequence 
of turnaround redshift in our model is a strict function of galaxy mass. 
The metallicities are well described by this model, including the lack of 
strong correlation of metallicity with redshift at the current level of obser- 
vational accuracy (( pettini et al., 1997 ); our Fig. ^3). Despite our deficit of 



low metallicity systems, we predict that there are numerous absorbers with 
metallicities in the range —2.0 < [Zn/H] < —0.8 waiting to be discovered 
with column densities 19.5 < log {Nhj cm~^) < 20.5. We also find that all 
our model absorbers have internal dust opacity r < 0.5 (with 50% below 
r = 0.1). This result can be interpreted in two very different ways. Taking 
our current model at face value, we predict that there are no significant 
selection effects in optical quasar surveys due to obscuration in foreground 



galaxies (see § 4.1). However, the current model does not consider interact- 
ing, starbursting systems which are expected to have much higher opacity, 
and which we will examine in a future paper ((Paker et al., 2001)). 



[INSERT FIGURE 5 HERE] [INSERT FIGURE 6 HERE] 



6.2.2 Global Properties of the Universe 

The evolution of the global properties of the Universe (as previously caL 
culated e.g. by ( [Wolfe et al., 1995 ; Storrie-Lombardi et al., 1996b[ ; Steidel 



3t al., 19991 ) ) emerge 'for free' from our model. The star formation history 



is shown in Fig. [8[, with data points derived from UV observations (as sum- 
marised by ( pteidel et al., 1999 ), without attempted dust corrections) shown 



for comparison with our independently derived smooth prediction. Our sim- 
ple model is a remarkably good representation of all the data at redshifts 
z < 4, but the real Universe appears to form more stars more quickly at 
the earliest observable times. Therefore, we would expect our high redshift 
{z > 4) metallicities to underpredict the observations, which indeed they 
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do, as shown in Fig. |^. 

[INSERT FIGURE 7 HERE] 

Our detailed prediction of the star formation history depends upon the 
adopted Holmberg relation and star formation parameter values. The obser- 
vations are at optical and UV wavelengths, which are particularly sensitive 
to steady, normal star formation or star formation bursts with low dust 
obscuration, and are strongly biased against heavily dust-embedded star 
formation, as seen in interacting and starbursting systems. Such embed- 
ded stellar energy generation is certainly comparable to that detectable at 
UV/optical wavelengths (( Puget and et al., 1996| ; Hauser and et al., 1998| ), 



and may be particularly important at higher redshifts {z > 1). Our mod- 
els do not include any effects of gas inflows, such as those triggered by 
interactions and mergers of galaxies, which appear to be the driving force 
behind dust embedded starburst activity. Therefore, it is particularly inter- 
esting that the global star formation history of our model (with steady, non- 
bursting star formation) is broadly similar, for redshifts z < 3, to the best 
UV/optical observations currently available. We compare our calculations 
to the UV/optical data without the uncertain observational dust reddening 
corrections, both because these corrections are so uncertain, and because 
they are likely to be small in UV-bright star formation regions. Galaxy 
evolution calculations which account for interactions, mergers and gas flows 
are required to study to true, total star formation history of the Universe. 

[INSERT FIGURE 8 HERE] 

In Fig. ^, we summarise the evolution of the baryons in our model clouds, 
as distributed between gas, stars, metals and dust (Fig. y). Our prediction 
of the total gas contribution in bound systems (dashed line) is in reasonable 
accord with the data (points from ( ^torrie-Lombardi et al., 1996b| )), which 
are based upon observations of HI in DLAs. Again, the data are consistent 
with more star formation at z > 4. 

[INSERT FIGURE 9 HERE] 



6.3 Detailed fits to individual DLAs 

We now move from broad-brush, statistical comparisons (§Q^) to detailed 
study of individual DLA features. We can only model systems in which the 
absorption redshift, Zabs^ the HI column density, Nhi and the metallicity 
([Z/H], as measured by [Zn/H], or [Fe/H] corrected for <[Zn/Fe]>£)LA) are 
known. We then take a slice in redshift space though our simulated uni- 
verse at Zahs (the most precisely determined observable) and extract all the 
galaxies with Nhi and [Z/H] which are consistent with the observations. 
This gives sets of model galaxy parameters, including turnaround redshifts. 
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masses, and impact parameters to the QSO line-of-sight. The results of our 
galaxy chemical evolution models for each observed DLA yield detailed star 
formation histories, which can then be used with a stellar spectral synthesis 
code ( (iBruzual and Chariot, 200C| )) to predict observable galaxy properties 
at the DLA absorption redshift, and indeed, any other epoch. 

We can now select a well-defined sample of DLAs for further observa- 
tional followup. Of the 63 DLAs with complete {zabsi ^Hi, [Z/H]) data, we 
could simultaneously obtain good fits to the observed redshift, column den- 
sity and metallicity of the DLA galaxy in 18 (~ 30 it 10%) cases (i.e. within 
the quoted (1 a) errors). It is likely that many of the other (unfitted) DLAs 
have undergone significant bulk gas fiows, possibly involving merger-driven 
star formation events. For those DLAs which we can fit well, the galaxy 
masses are M ~ 3 x 10^ - 3 x 10^'^ Mq, the galaxies are at impact parame- 
ters of 6 < 2kpc, and will have apparent magnitudes of 17 < i^ < 28. We 
therefore predict that a large fraction of DLAs are caused by sub-L* galax- 
ies sitting almost directly in front of the QSO. This prediction implies that 
the ~ L* galaxy associated with the absorber Q2233+1310 {zabs = 3.151) 
by Steidel and others (( [Steidel et al., 1995 ; Djorgovski et al., 1996[ )), about 



15h~^ kpc from the quasar sight line (in the same cosmology as our individ- 
ual model fits) is either atypical, or is not the true absorber. 

For each DLA we can fit well, we have selected up to three candidate 
model galaxies which bracket the plausible mass distribution, and evolved 
these from formation to their absorption redshift, using stellar synthesis to 
estimate key observables such as the galaxy K-band apparent magnitude (in 
the appropriate cosmology). The predicted galaxy properties for each DLA 
are summarised in Table IT], with typical uncertainties being ±0.5 kpc in the 
impact parameter and ±0.5 in the absolute and apparent magnitudes. 

[INSERT TABLE 1 HERE] 

As would be expected from our galaxy formation prescription, galaxy 
luminosity is broadly anticorrelated with (absorption) redshift. We predict 
that some DLA galaxies will be too faint and too close to the quasar to 
be detectable even with near-IR detectors on 8m telescopes. But there 
is a small, accessible population with predicted K < 22, and b > 2 kpc at 
redshifts z < 2.5. The optical-infrared colours are generally moderate (B— K 
< 5.5), but up to a third of these DLAs arise in rather red (6 < B— K < 
8.5) objects, typically at redshifts 2 < Zabs < 3. 

These results demonstrate the importance of detailed galaxy chemical 
evolution modelling in the design of large telescope programmes to observe 
DLAs and galaxies. For example, ( Aragon-Salamanca et al., 1996| ) (AEO) 



use a 4m-class telescope to image 10 DLA quasars at K, and found two 
candidate DLA galaxies, towards Q0841 ± 129 (which has 3 DLA systems) 
and Q1215 ± 333, with i^ -^ 20 and b ~ 10 kpc; they set magnitude limits 
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of i^ ~ 20.3-20.7 for the others. We predict that the true absorbing galaxy 
towards Q1215 + 333 {zabs = 1-99) has a mass Mb ~ 5 x lO^M© , with 
X ~ 22 and impact parameter of only 6 ~ 1 kpc. The closest DLA towards 
Q0841 + 129 is at Zahs = 2.37, and we predict an absorbing galaxy with mass 
Mf, ~ 5 X 10^ M0 , with K ~ 21.5 and impact parameter of only h ~ Ikpc 
(the Zabs = 2.48 system will be very similar). The other DLAs in the AEO 
sample are also probably sub-L* galaxies at small impact parameters, well 
beyond the apparent magnitude and angular resolution limits of AEO. We 
therefore do not believe that AEO have (nor could have) detected the true 
DLA absorbers in the DLAs they studied. Without spectroscopic galaxy 
redshifts, we can only speculate that the AEO galaxies may be in the same 
cluster as the true DLA absorbing galaxies. 

6.4 Local Group dwarf galaxies 

We can do 'backyard cosmology' by comparing the properties of the dwarf 
galaxies in the Local Group with our models. If we believe that dwarf 
galaxies give rise to a significant fraction of the observed DLAs, then we must 
explain why many dwarfs in the Local Group are too gas-poor to produce 
the column densities in DLAs. This may be an environmental effect. The 
Local Group dwarf galaxies clearly separate into two classes ( ([Mateo, 1998]) ): 
gas-poor dwarfs, those within ~ 250 kpc of the Milky Way or Andromeda; 



and gas-rich dwarfs, which lie further out, as shown in Fig |10|a. Dwarfs near 
giants have probably had their gas stripped by tidal interactions. 

The stellar metallicities in the Local Group dwarf galaxies are the mean 
values measured for old and intermediate-age stars (Fig |lO|b), and thus re- 
flect the gas-phase metallicities 1-10 Gyr ago. The observed values, —2.0 < 
[Fe/H] < —0.5, are similar to the properties of DLAs at redshifts z ~ 2, 
which in our low density cosmology corresponds to a look-back time ~ 6 
Gyr. There may be a hint that the lowest mass dwarf galaxies (< 10^ Mq ) 
have the lowest metallicities, although the circum-Galactic dwarfs dominate 
the statistics. Low mass galaxies may be subject to the suppression of star 
formation and hence chemical evolution due to supernova winds, leaving 
gas-poor, low metallicity systems. Alternatively, low mass galaxies may 
simply have very low star formation rates, giving a gas-rich, low metallicity 
system. The gas fraction of the majority of Local Group dwarf galaxies are 
comparable to those inferred for DLA absorbers from our models. 

[INSERT FIGURE 10 HERE] 

We therefore conclude that there are dwarf galaxies in the bulk of the 
volume of the Local Group which have gas fractions and metallicities com- 
patible with those of DLAs. Further investigation would be desirable, but 
we do not attempt the necessary more detailed modelling of Local Group 
dwarfs here. 
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6.5 Time delayed elements 

The main sources of heavy elements are thought to be Type II supernovae, 
which enrich the interstellar medium 'promptly' on timescales of ~ 5 x 10^ 
yr, and Type la supernovae, which are 'delayed' by up to ~ 1 Gyr. The 
elements predominantly produced in massive stars are thought to include 
Si, O, Mn and Zn. It is thought that Fe is partly produced in SNII, but 
that SNIa may be at least twice as efficient Fe producers. We have therefore 
modelled the evolution of element ratios by calculating a delayed metallicity 
as well as the standard prompt metallicity which we compare with [Zn/H] 
measurements throughout this paper. The evolution of 'Promptium' (which 
can masquerade as e.g. O, S, Si) with respect to an element, 'Delaydium', 
which is one third prompt and two thirds delayed {e.g. Fe) can be calculated 
for different relative time delays, such as 0.25, 1 and 3 Gyr. 

[INSERT FIGURE 11 HERE] 

At early times (high redshift, low metallicity) in the model systems, only 
the prompt portion of iron has been produced, and so [Pr/De] = log ( |) = 
0.47. After r = 1 Gyr, the delayed portion of iron starts to build up and 
[Pr/De] falls towards the final cosmic value, [Pr/De] = 0. To reach higher 
values such as [Pr/De] ~ 1 would require that Delaydium only has e.g. ~ 10% 
prompt component. 

In Fig. |l^, we compare the model predictions for the build-up of iron 
(Delaydium) relative to oxygen, sulphur and silicon (Promptium) for a time 
delay r = 1 Gyr in the low density cosmology with the best currently avail- 
able data, as a function of (a) absorption redshift and (b) metallicity as 
measured by iron (Delaydium). These data are clearly subject to the usual 
uncertainties about depletions, and the putative prompt production path- 
ways. Ironically, the most theoretically useful prompt element is oxygen, 
for which the observations are still subject to particularly significant ambi- 
guities (( Prochaska and Wolfe, 20001) ). At this stage, all we can say is that 



our synthetic DLA survey predictions are in reasonable agreement with the 
available data. More informative investigations of element ratios must await 
abundance measurements with an accuracy <7[z/Hl ~ 0.1, which are starting 
to become available (( prochaska and Wolfe, 2000| )). 

7 Limitations of the Model 

We have made no attempt in this paper to model the role of mergers in 
the formation and evolution of galaxies. However, it is clear from infrared 
surveys (which are efficient at finding interacting galaxies ( (|Clements et al.~ 



1996| ))), and from theoretical work ( ([Barnes and Hernquist, 1996| )), that 



mergers are a fundamental galaxy formation process. Why is it therefore 
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that our model galaxies do so well in reproducing the observational data for 
DLAs? As a step towards answering this question, we outline a scenario in 
which mergers do play a role. 

Suppose that the first generation of bound objects in the Universe were 
small sub-galactic units of mass ~ 10^ - 1{)^ Mq which somehow produce 
stars. These low-mass star forming units are likely to lose their gas in 
supernova-driven winds after only a few generations of star formation. If 
we allow some of these gas-poor, low-mass stellar systems to merge, they 
would form larger bound systems, embedded in halos of lightly processed 
ejected gas. A fraction of this lightly enriched halo gas could settle into a 
disk around the merged stellar systems, and more familiar population I and 
II star formation could proceed. 

The statistics of the merging systems and the subsequent synthesis of 
giant galaxies can be described by the Press-Schechter model of galactic 
structure formation ( ([Press and Schechter, 1974 )). The interesting point 



about the scenario outlined above, is that only a small amount of star for- 
mation takes place prior to the formation of the disk galaxy. In our simple 
model, very little star formation takes place during the collapse phase. This 
suggests that the present day state of a giant galaxy as predicted by this type 
of model may be insensitive to the details of its formation mechanism ( ([Katz 



and Gunn, 199]] )), and hence, our preliminary model can still appropriately 
describe the data. 



8 Conclusions and Next Steps 

We have constructed a simple, robust model of galaxy chemical evolution, 
which we quite successfully applied to the study of galaxies detected as 
damped Lya systems in the spectra of quasars. We have modelled the 
evolution of neutral hydrogen column density, metallicity and hence, dust 
opacity, in clouds of initially primordial gas from the point where they begin 
to collapse under self gravity, up to the present. We find that the outer re- 
gions of model giant galaxies, and model sub-L* galaxies are able to produce 
absorbers with column density and metallicity typical of DLAs, and remain 
optically thin. We have studied the specific properties of individual DLAs 
and find that the models of the underlying galaxies are in good agreement 
with the limited imaging observations which currently exist, making sense 



of the relative lack of success in detecting DLAs {e.g. ( Ar agon- Salamanca 



|et al., 1996P ). Contrary to some expectations (( [Wolfe et al., 19951 ) ), we 
predict that a significant proportion of the DLA systems at moderate red- 
shifts {z ~ 2.5) arise in sub-L* galaxies. We have therefore initiated a new 
programme of infrared imaging using the Gemini North telescope, which 
is capable of detecting these fainter galaxies at smaller impact parameters 
around DLA quasars. 
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Our 'monolithic' models provide a reasonable account of the chemical 
evolution of the Universe at redshifts below z ~ 3. However, in attempt- 
ing to use the Holmberg relation to constrain the cosmic history of galaxy 
formation, we have exposed the limitations of our model assumption con- 
cerning the unimportance of galaxy mergers, particularly in promoting star 
formation at high redshifts (2 > 3). Our next step is therefore to a slightly 
more sophisticated galaxy formation scenario, including mergers, to search 
for observable effects, including changes in predicted DLA properties. 
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Figure Captions 

Figure 1 The evolution of the quasar optical luminosity function (QOLF), 
and the effects of obscuration, (a) A parametric fit to the QOLF at red- 
shift z = 2; (b) A pure luminosity parameterisation of the evolution of the 
characteristic quasar luminosity (M^); (c) The apparent magnitude of the 
'knee' in the QOLF, accounting for PLE and cosmological dimming (error 
bar indicates approximate overall uncertainties); (d) The foreground optical 
depth required to dim an L*^ quasar below a survey apparent magnitude 
limit friBn^ = 19 (and uncertainties). 

Figure 2 Typical dust opacity predicted from the gas phase metallicity 
and the neutral gas column density. The dashed lines correspond to typical 
dust optical depths r^ = 1, 0.1, 0.01. The solid box shows the region occupied 
by normal galaxies (with a possible extension into the dotted box) . The data 
points at [-1.8, 21.8] and [-1.6, 21.8] are I Zw 18 and the LMC. The remaining 
data points are all known DLAs for which metallicity measurements are 
available ([Fe/H] measurements have been corrected for <[Zn/Fe]>£)Lyi (see 

Figure 3 The model redshift selection function, which is a fit to the com- 
bined redshift survey path lengths, X(z), of the existing large, well defined 
DLA surveys (see text). 

Figure 4 Comparison of DLA observations (points) with the predictions 
of our model (contours) using Monte Carlo simulations in Einstein-de Sit- 
ter cosmology. The dashed lines indicate the approximate selection effects 
which affect the data (derived from DLA surveys, including the metallicity 
measurements of Pettini and collaborators; see text). 

Figure 5 Cumulative distributions of observed and synthetic DLAs in 
(a) redshift, (b) column density and (c) metallicity in the low density uni- 
verse. Identical selection effects, empirically derived from DLA surveys (in- 
cluding the metallicity measurements of Pettini and collaborators; see text), 
have been applied to the model and the data. 

Figure 6 Comparison of DLA observations (points) with the parameter 
space density of systems predicted by our model (contours) using Monte 
Carlo simulations in low density cosmology. Identical selection effects, em- 
pirically derived from DLA surveys (including the metallicity measurements 
of Pettini and collaborators; see text), have been applied to the model and 
the data. 

Figure 7 Comparison of DLA metallicities as a function of redshift from 
observations (points, [Zn/H] and [Fe/H]) and our model (contours), as in 
Figure 6 hut with no selection effects applied. 

Figure 8 The star formation history of the Universe. The data points 
are for UV observations, taken from Steidel 1999 and not corrected for ex- 
tinction. Our model prediction is shown as a (smoothed) curve. The curve 
has not been fitted to these data. 
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Table 1: Predicted properties of the absorbing galaxy from successful fits 
to individual DLAs. (Assumed cosmology: Q^m = 0.3, il^ = 0.7, Q.\ = 0.0 
('low density')). 



DLA 


Absorption 


Impact 


Mass 


Mass 








name 


redshift 


parameter 


(baryons) 


(stars) 










^ahs 


h (/i-^kpc) 


log(M0 ) 


log(M0 ) 


Mm 


niK 


B-K 


0100+13000 


2.309 


0.8 


10.5 


8.0 


-22.2 


19.4 


5.4 


0112+03000 


2.423 


0.0 


7.9 


6.0 


-15.3 


27.6 


8.5 


0149+33500 


1.613 


1.4 


10.0 


7.8 


-20.6 


21.1 


5.5 


0454+03600 


0.860 


0.9 


10.0 


8.3 


-19.0 


22.3 


5.5 


0528-2500A 


2.140 


0.8 


9.7 


7.6 


-20.0 


21.8 


5.6 


0841+1256A 


2.375 


0.6 


9.6 


7.4 


-19.9 


22.0 


5.7 


0841+1256B 


2.476 


0.9 


9.8 


7.5 


-20.6 


21.2 


5.6 


1151+06800 


1.774 


1.7 


11.2 


8.8 


-23.5 


17.8 


5.4 


1215+33300 


1.999 


0.6 


9.7 


7.6 


-19.6 


22.1 


5.6 


1223+17530 


2.466 


0.4 


10.3 


7.9 


-22.3 


19.4 


5.3 


1229-02100 


0.395 


0.4 


9.3 


7.9 


-17.0 


22.8 


5.5 


1247+26700 


1.223 


0.4 


8.4 


6.8 


-15.7 


26.4 


6.0 


1328+30700 


0.692 


2.8 


11.5 


9.7 


-23.3 


17.4 


5.5 


1331+17040 


2.800 


0.7 


10.2 


8.1 


-20.8 


21.5 


7.1 


2206-1990B 


2.076 


4.0 


10.9 


8.5 


-23.3 


18.1 


5.2 


2344+12400 


2.538 


0.1 


7.8 


5.9 


-15.1 


27.1 


6.2 


2348-14440 


2.279 


1.1 


9.9 


7.7 


-20.7 


21.1 


5.6 


2359-0220B 


2.152 


0.5 


8.9 


7.0 


-17.6 


24.3 


5.8 
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Figure 9 The global properties of our model Universe. The fraction of 
the closure mass density supplied by baryons in galaxies, ^b^„ (solid line); 
gas in galaxies, ^gas^„ (dashed) ; in stars, Q.stars (dotted); in metals, ^metals 
(dash-dot), in dust, i^dust (dash-dot-dot-dot). The contribution from to- 
tal mass in galaxies, ^toti„ = lO^fein- The contribution from mass not yet 
formed into galaxies, ^tot,out = ^M — ^toUni where the overall mass con- 
tribution Qm = 0.3 in our low density cosmology. The data points are the 
inferred ^gas in DLAs (([Storrie-Lombardi et al., 1996a; Rao and Turnshek, 



2000| ; IStorrie-Lombardi and Wolfe, 200q )). 



Figure 10 Local Group dwarf galaxies as a function of distance from 
the Milky Way (a) HI gas fraction (b) Stellar metallicity [Fe/H] for dwarf 
galaxies with mass < 10^ Mq (open circles) and > 10^ Mq (crosses). 

Figure 11 Ratios of prompt element (Promptium, Pr) to delayed el- 
ement abundances (Delaydium, De) for simulated DLA survey (points as 
dots), with data for [0/Fe] (circles), [S/Fe] ('S' shapes) and [Si/Fe] (tri- 
angles), with errors) as a function of a) absorption redshift b) metallicity 
([De/H] or [Fe/H]). 
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